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Early indications by H.E.S.S. and the subsequent detection of blazar 3C279 by MAGIC 
show that the Universe is more transparent to very-high-energy gamma rays than previously 
thought. We demonstrate that this circumstance can be reconciled with standard blazar emis- 
sion models provided that photon oscillations into a very light Axion-Like Particle occur in 
extragalactic magnetic fields. A quantitative estimate of this effect indeed explains the ob- 
served spectrum of 3C279. Our prediction can be tested by the satellite-borne Fermi/LAT 
detector as well as by the ground-based Imaging Atmospheric Cherenkov Telescopes H.E.S.S., 
MAGIC, CANGAROO III, VERITAS and by the Extensive Air Shower arrays ARGO-YBJ 
and MILAGRO. 



1 Introduction 

A characteristic feature of the very-high-energy (VHE) band is that the horizon of the ob- 
servable Universe rapidly shrinks above 100 GeV as the energy further increases. This is due 
to the fact that photons from distant sources scatter off background radiation permeating the 
Universe, thereby disappearing into electron-positron pairs ^4 The corresponding cross sec- 
tion 17(77 —> e + e~) peaks where the VHE photon energy E and the background photon en- 
ergy e are related by e ~ (500 GeV/E) eV. We recall that Imaging Atmospheric Cherenkov 
Telescopes (IACTs) probe the energy interval 100 GeV — 100 TeV. Consequently, observa- 
tions performed by the IACTs are affected by an opacity dominated by the interaction of 
the beam photon with ultraviolet /optical/infrared diffuse background photons (frequency band 
1.2 • 10 3 GHz — 1.2 • 10 6 GHz, corresponding to the wavelength range 0.25 fim — 250 yum), usually 
called Extragalactic Background Light (EBL) and produced by galaxies during the whole his- 
tory of the Universe. Neglecting evolutionary effects for simplicity, photon propagation is then 
controlled by the photon mean free path X^(E) for 77 — > e + e~, and so the observed photon 
spectrum ^ Q ^ S (E, D) is related to the emitted one $> em (E) by 



$ obs (£,D) = e- D / A ^) <J> em (£) 



(1) 



Within the considered energy range, X^(E) decreases like a power law from the Hubble 
radius 4.3 Gpc around 100 GeV to nearly IMpc around 100 TeV^. Thus, Eq. ((TJ implies that 
the observed flux is exponentially suppressed both at high energies and at large distances, so 
that sufficiently far-away sources become hardly visible in the VHE range and their observed 
spectrum should anyway be much steeper than the emitted one. 

Yet, the behaviour predicted by Eq. ([TJ has not been detected by observations. A first 
indication in this direction was reported by the H.E.S.S. collaboration in connection with the 
discovery of the two blazars H2356-309 (z = 0.165) and 1ES1101-232 (z = 0.186) at E ~ 1 TeV^l. 
Stronger evidence comes from the observation of blazar 3C279 (z = 0.536) at E ~ 0.5 TeV by 
the MAGIC collaboration SI. In particular, the signal from 3C279 collected by MAGIC in the 
region E < 220 GeV has more or less the same statistical significance as the one in the range 
220 GeV < E < 600 GeV (6.1a in the former case, 5.1a in the latter) |. 

A possible way out of this difficulty involves the modification of the standard Synchro-Self- 
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Compton (SSC) emission mechanism. One option invokes strong relativistic shocks^. Another 
is based on photon absorption inside the blazar^l. While successful at substantially hardening 
the emission spectrum, these attempts fail to explain why only for the most distant blazars does 
such a drastic departure from the SSC emission spectrum show up. 

Our proposal - usually referred to as the DARMA scenario - is quite different^. Implicit 
in all previous considerations is the hypothesis that photons propagate in the standard way 
throughout cosmological distances. We suppose instead that photons can oscillate into a new 
very light spin-zero particle - named Axion-Like Parlicle (ALP) - and vice-versa in the presence 
of cosmic magnetic fields, whose existence has definitely been proved by AUGER observations^. 
Once ALPs are produced close enough to the source, they travel unimpeded throughout the 
Universe and can convert back to photons before reaching the Earth. Since ALPs do not undergo 
EBL absorption, the effective photon mean free path X^^g(E) gets increased so that the observed 
photons travel a distance in excess of X-y(E). Correspondingly, Eq. ([1]) becomes 

$ ohB {E,D) = e -D/^ME) $ cm( £) ; (2) 

which shows that even a small increase of Xy >e g(E) gives rise to a large enhancement of the 
observed flux. It turns out that the DARMA mechanism makes X~, te g(E) shallower than X~(E) 
although it remains a decreasing function of E. So, the resulting observed spectrum is much 
harder than the one predicted by Eq. (pQ) , thereby ensuring agreement with observations even for 
a standard SSC emission spectrum. As a bonus, we get a natural explanation for the fact that 
only the most distant blazars would demand $> em (E) to substantially depart from the emission 
spectrum predicted by the SSC mechanism. 

We proceed to review the main features of our proposal as well as its application to blazar 
3C279. 

2 DARMA scenario 

Both phenomenological and conceptual arguments entail that the Standard Model (SM) of 
particle physics should be viewed as the low-energy manifestation of some more fundamental 
and richer theory of all elementary-particle interactions including gravity. Therefore, the SM 
lagrangian is expected to be modified by small terms describing interactions among known and 
new particles. Many extensions of the SM which have attracted considerable interest in the last 
few years indeed predict the existence of ALPs. They are spin-zero light bosons defined by the 



"See ref.^ for a different view. 



low-energy effective lagrangian 



£alp = ^"a^a-^^^—r^a, (3) 
where F^ v is the electromagnetic field strength, Fn V is its dual, a denotes the ALP field and m 

— 1/2 

stands for the ALP mass. According to the above view, it is assumed M » G F ~ 250 GeV. 
On the other hand, it is supposed that m <C G F 1 ^ 2 ~ 250 GeV. The standard AxionQ^I i s th e 
archetype of ALPs and is characterized by a specific relation between M and m, while in the 
case of generic ALPs M and m are to be regarded as independent. So, the peculiar feature of 
ALPs is the trilinear 7-7-a vertex described by the last term in Calp, whereby one ALP couples 
to two photons. 

Owing to such a vertex, ALPs can be emitted by astronomical ob jects of various kinds, and 
this fact yields strong bounds: M > 0.86 • 10 10 GeV for m < 0.02 eV^ and M > 10 11 GeV for 
Moreover, the same 7-7-a vertex produces an off-diagonal element in the mass 
matrix for the photon-ALP system in the presence of an external magnetic field B. Therefore, 
the interaction eigenstates differ from the propagation eigenstates and photon-ALP oscillations 
show upQ21. 

We imagine that a sizeable fraction of photons emitted by a blazar soon convert into ALPs. 
They propagate unaffected by the EBL and we suppose that before reaching the Earth a substan- 
tial fraction of ALPs is back converted into photons. We further assume that this photon-ALP 
oscillation process is triggered by cosmic magnetic fields (CMFs), whose existence has been 
demonstrated very recently by AUGER observations^. Owing to the notorious lack of informa- 
tion about their morphology, one usually supposes that CMFs have a domain-like structure^! 
That is, B ought to be constant over a domain of size L^ ora equal to its coherence length, with 
B randomly changing its direction from one domain to another but keeping approximately the 
same strength. As explained elsewhere^', it looks plausible to assume the coherence length in 
the range 1 — lOMpc. Correspondingly, the inferred strength lies in the range 0.3 — 1.0 nG. 



3 Predicted energy spectrum 

Our ultimate goal consists in the evaluation of the probability P~->y(E, D) that a photon remains 
a photon after propagation from the source to us when allowance is made for photon-ALP 
oscillations as well as for photon absorption from the EBL. As a consequence, Eq. ([2]) gets 
replaced by 

$ obs (E, D) = P 7 ^ 7 (E, D) $ em (E) . (4) 

Our procedure is as follows. We first solve exactly the beam propagation equation arising from 
£alp over a single domain, assuming that the EBL is described by the "best-fit model" of 
Kneiske et a/.^l Starting with an unpolarized photon beam, we next propagate it by iterating 
the single-domain solution as many times as the number of domains crossed by the beam, taking 
each time a random value for the angle between B and a fixed overall fiducial direction. We 
repeat such a procedure 10' 000 times and finally we average over all these realizations of the 
propagation process. 

We find that about 13% of the photons arrive to the Earth for E = 500 GeV, representing 
an enhancement by a factor of about 20 with respect to the expected flux without DARMA 
mechanism (the comparison is made with the above "best-fit model"). The same calculation 
gives a fraction of 76% for E = 100 GeV (to be compared to 67% without DARMA mechanism) 
and a fraction of 3.4% for E = 1 TeV (to be compared to 0.0045% without DARMA mechanism). 
The resulting spectrum is exhibited in Fig. 1. The solid line represents the prediction of the 
DARMA scenario for B ~ 1 nG and L^ oxa ~ 1 Mpc and the gray band is the envelope of the 



results obtained by independently varying B and Ldom within a factor of 10 about such values. 
These conclusions hold for m < 10~ 10 eV and we have taken for definiteness M ~ 4 • 10 11 GeV 
but we have cheked that practically nothing changes for 10 11 GeV < M < 10 13 GeV. 

Our predictions can be tested by the satellite-borne Fermi/LAT detector as well as by the 
ground-based IACTs H.E.S.S., MAGIC, CANGAROO III, VERITAS and by the Extensive Air 
Shower arrays ARGO-YBJ and MILAGRO. 




Figure 1: The two lowest lines give the fraction of photons surviving from 3C279 without the DARMA mechanism 
within the "best-fit model" of EBL (dashed line) and for the minimum EBL density compatible with cosmology 
(dashed-dotted line), which are discussed by Kneiske et al. (see reference in the text). The solid line represents 

the prediction of the DARMA mechanism. 
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